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Abstract 

We have studied theoretically the electron spin relaxation in semiconductor 

quantum dots via interaction with nuclear spins. The relaxation is shown 

to be determined by three processes: (i) - the precession of the electron 

spin in the hyperfine field of the frozen fluctuation of the nuclear spins; (ii) 

- the precession of the nuclear spins in the hyperfine field of the electron; 

and (iii) - the precession of the nuclear spin in the dipole field of its nuclear 

neighbors. In external magnetic fields the relaxation of electron spins directed 

along the magnetic field is suppressed. Electron spins directed transverse to 

the magnetic field relax completely in a time on the order of the precession 

period of its spin in the field of the frozen fluctuation of the nuclear spins. 

Comparison with experiment shows that the hyperfine interaction with nuclei 

may be the dominant mechanism of electron spin relaxation in quantum dots. 
PACS: 72.25.Rb, 78.67.Hc, 78.67.Bf 
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I. INTRODUCTION 



The long electron spin dephasing time (hundred nanoseconds) recently reported sug- 
gests using the spin of an electron localized in a quantum dot (QD) as the realization of a 
quantum bit, and electron doped quantum dots as the basic material for implementing a 
solid state based quantum computer 0. Room temperature operation, which is usually a 
primary requirement in microelectronics, is of secondary importance for quantum computers 
relative to that of finding a material for their physical realization. Rather, operating now at 
low temperature, where the localized electron has a long spin coherence time, is an essential 
condition for quantum computing and spin memory storage. The dominant electron spin re- 
laxation mechanism in bulk is connected with the spin-orbit interaction of carriers (see Ref. 
||), but is strongly suppressed for localized carriers Q. Spin relaxation due to the electron 
hole exchange interaction plays an important role during the time of nonequilibrium carrier 
relaxation but does not effect localized electron spin relaxation afterword. Dephasing 

of the electron in QD ground states via two phonon real transitions to, or virtual transitions 
through, excited states (Urbach or Raman two phonon processes ||) is also suppressed at 
low temperatures || . As a result electron spin relaxation via interaction with nuclei becomes 
the dominant relaxation mechanism for localized electrons at low temperatures. 

The interaction of localized electrons with nuclei was studied early on for electrons lo- 
calized at donors in bulk (see for example fllOlJTTll ). There, the electron interacts with a 
large number of nuclei and feels the hyperfine magnetic field of the nuclei located in the 
region where electron is localized; this is also true for QD localized electrons. However, the 
correlation time of the electron-nucleus interaction of donor localized electrons is limited by 
the time of shallow donor ionization, and tunnel jumps between the donors, r c . Usually this 
time is much shorter than the period of the electron spin precession, to J , in the hyperfine 
field of the nuclei. For quantum computation and spin storage, however, it is important 
to have a large value of r c , and we are interested in the opposite limit: 0JfT c ^> 1. This is 
realized in QDs at sufficiently low temperatures. 
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In this paper we consider the electron spin relaxation via its interaction with the spins of 
the nuclei in QDs in the absence and the presence of an external magnetic field. The process 
is facilitated by the disparity of the characteristic time scales of the three processes that 
determine the relaxation: the period of the electron precession in the frozen fluctuation of 
the hyperfine field of the nuclei, the period of the nuclear spin precession in the hyperfine field 
of the electron, and the nuclear spin relaxation time in the dipole-dipole field of its nuclear 
neighbors. Estimates of these time scales can be made for the case of GaAs, whose hyperfine 
constants are well k own. For QDs containing 10 5 nuclei they are found to be: ~ 1 ns, ~ 1 /is 
and ~ 100 /is, respectively. Therefore, as a first step, we can describe the electron spin 
relaxation as a precession in the quasi-stationary frozen fluctuation of the hyperfine field 
of the nuclear spins. We can than examine the additional electron spin relaxation arising 
from the slowly varying nuclear hyperfine fields. The nuclear dipole-dipole interaction does 
not conserve the total nuclear spin; the third time scale provides a natural limit on the 
coherence of the electron-nuclei spin system. However many other relaxation mechanisms 
are important at this long time scale, therefore and we will not consider the effect of the 
nuclear dipole-dipole interactions on the electron spin relaxation. A spin dephasing time 
(T 2 * ) calculated for an ensemble of QDs is in a good agreement with available experimental 
data. 

The paper is organized as follows: The hyperfine interaction of an electron with the nuclei 
and the hyperfine nuclear magnetic fields acting on the electron are considered in section 
II. The electron spin relaxation times in the absence and in the presence of an external 
magnetic field are calculated in sections III and IV respectively. In section V we compare 
our theoretical results with experimental data. 
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II. HYPERFINE INTERACTION OF A LOCALIZED ELECTRON WITH 
NUCLEI. THE FROZEN FLUCTUATION OF THE NUCLEAR HYPERFINE 

FIELD IN QUANTUM DOTS. 



The electron spin relaxation due to the nuclear spins is determined by their hyperfine 



Fermi contact interaction. The Hamiltonian of this interaction can be written ||12 

16n ^ fMj 



H, 



cont 



^Y,f(S-I 3 )S(r-Rj) , (1) 
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where //g is the Bohr magneton, S and r are the spin and position of the electron, fij, Ij 
and Rj are the magnetic moment, spin and position of the j-th nucleus, and the sum goes 
over all the nuclei in the lattice. For localized electrons the distance between their energy 
levels is much larger than the energy of the hyperfine interaction with the nuclei. As a result 
this spin-spin interaction can be described by the Hamiltonian: 

Hu = %E Ai WW\ 2 {H** + &* + H*v) . ( 2 ) 

Z 3 

which is obtained from first order perturbation theory. Here Vq is the volume of the unit 
cell, ip{Rj) is the electron envelope wave function at the j-th nucleus, I a and a a are the spin 
projections on the coordinate axes a = x,y,z, and A 7 = {lQ'nn B iij/'5P)\u c (Rj)\ 2 , where 
u c (Rj) is the electron Bloch function at the nucleus. In GaAs, the sum of A? over all the 
nuclei in the unit cell A = Y,j A? ~ 90 /xeV [0 . We can neglect the interaction of holes with 
the nuclei because the hole Bloch functions vanish at the nuclear positions. In addition, we 
neglect, for now, the nuclear dipole-dipole interactions, which do not conserve the total spin 
of the electron-nuclear system. They become important only at times longer than 10 -4 s. 

The effective nuclear hyperfine magnetic field, Bn, acting on a localized electron spin 
can be obtained from Eq.|^ taking the expectation of the Hamiltonian over the ensemble 
of nuclear wave functions. It is the sum of contributions from a large number of nuclei: 

B N = -^(£A>MR j )\*?) N (3) 
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where (...)n denotes a quantum mechanical average over the ensemble of nuclear wave func- 
tions and g e is the electron ^-factor. We should note that the precession frequency of the 
electron in the hyperfine field of all the nuclei is much greater then the precession frequency 
of a nucleus in the hyperfine field of the electron. That is, the electron sees a snapshot of 
the "frozen fluctuation" of the nuclear field. The magnitude and direction of this field are 
randomly distributed, and described by a Gaussian probability density distribution function: 



A! 



(4) 



where A# is the dispersion of the nuclear hyperfine field distribution, 



A| = ^((B iV ) 2 ) = ^P(P + l)(a J ) 2 , (5) 

where 

a, = (W^ e )^X^)| 2 (6) 

is the magnetic field of a single nuclear spin acting on the electron, and we assumed that 
the nuclear spin directions are independent of each other. All nuclei in GaAs have the same 
spin P — I — 3/2. Replacing the sum over unit cells by an integration we obtain: 



x ■ 2J(/+l)E J (^) 2 ^o 16/(7 + 1)£jM 2 
B 3 (fi B g e y V L 3N L (fi B g e y ' 1 } 

where the sum in this equation goes over only those nuclei in a unit cell, 

VL = (/ d z r^\r))~ l and N L = 8V L /v (8) 

is the number of nuclei in the volume, VL, that effectively determine the electron precession 
frequency. In GaAs, the sum J2j{^) 2 ~ 1-2 • 10~ 3 meV 2 . The volume, VL, is on the order of 
the volume of the electron localization. 



III. ELECTRON SPIN RELAXATION IN ZERO EXTERNAL MAGNETIC FIELD 
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A. Electron spin dephasing in the frozen fluctuations of the nuclear field 



Let us consider an ensemble of identical QDs in which we simultaneously (at time t — 0) 
create electrons all having the same spin orientation (They can be created, for example, 
by circularly polarized light). The nuclear spins in the QDs of this ensemble are randomly 
oriented, the nuclear hyperfine fields in the dots differ from one another and, therefore, have 
a different effect on the initial electron spin, Sq, in each dot. We will consider the time 
dependence of the ensemble average electron spin relaxation for times small relative to the 
period of the nuclear precession in the hyperfme field of the electron. Each electron spin 
will be moving in the frozen fluctuation of the nuclear hyperfine magnetic field, Bn, (see 
Eq.|3|) in its own QD. These fields, however, are randomly distributed among the dots of the 
ensemble. Therefore, even though each electron spin will precess in a coherent fashion in the 
frozen hyperfine field of its own dot, the ensemble average spin polarization will decrease. 

The equation of motion of the spin S in a fixed magnetic field, B, is given by: 



S(t) = (S ■ n)n + {S — (S ■ n)n} cosut 
+ [{So — (So ■ n)n} x n] smu>t , 



(9) 



where So is the initial spin, n = B/B is a unit vector in the direction of the magnetic 
field, and u = fiBdeB/h is the Larmor frequency of the electron precession in this field. 
The equation also describes the coherent electron spin precession in a single QD due to 
the magnetic field, B^, of the frozen fluctuation of the nuclei (n = B^/B^, and uo = 
figg e B^/h). Averaging Eq.(|9|) over the magnetic field distribution of Eq.(|), we obtain the 
time dependence of the ensemble averaged electron spin polarization: 

t ^ 



(S(t)) = y{l + 2 



exp 



(10) 



The same time dependence describes the electron spin polarization of a single quantum dot 
averaged over a large number of measurements. Here 

n 



ft. 



3A T , 



(11) 



is the ensemble dephasing time which arises from the random electron precession frequencies 
in the randomly distributed frozen fluctuation of the nuclear hyperfme field in the dots. This 
time is on the order of 1 ns for GaAs quantum dots with 10 5 nuclei. The spin dephasing 
time is proportional to \/Vl- One can see that the average electron polarization relaxes to 
10% of its original value after a time equal to the dephasing time and then increases to a 
steady state value of 33% of its initial polarization. 

B. Electron spin dephasing as a result of variations of the nuclear field direction 

A localized electron interacts with a large number of nuclei, Nl 3> 1- The interaction 
of the electron spin with a single nucleus is ^/Wl times weaker than its interaction with the 
effective magnetic field of the frozen nuclear fluctuation. Changing the direction of a single 
nuclear spin only weakly perturbs the electron spin motion. The precession of the electron in 
the macroscopic fluctuation of the nuclear spins is s/Nl times faster than the precession of a 
nucleus in the hyperfme field of an electron, i.e., the nuclear precession period T N ~ T^N L . 
Since an electron in a QD precesses so rapidly around the nuclear magnetic field Bn, the 
nuclei only see the long time average of the hyperfme field of the electron, which is directed 
along Bn. Note, that components of the nuclear field perpendicular to B^ cancel each other 
out. Now, however, each nucleus, in turn, precesses about this direction with a different 
precession rate that is proportional to the square of the electron wave function at their 
respective nuclear positions. This variation in the precession rates leads to a non-vanishing 
slow time varying change in the frozen fluctuation of the hyperfme magnetic field of the 
nuclei, AB N ~ J2j (^IVK-^j)! 2 ) 2 x Bn)n, that is perpendicular to B N . These random 
changes in the nuclear magnetic field result in an additional relaxation of the electron spin 
polarization. 

Equation ( |T0"D describes the time dependence of the electron spin relaxation in the frozen 
fluctuation of the nuclear hyperfine fields for times much less than T/v- To include the effect 
of the time dependent changes in the nuclear fields, we need to examine the ensemble average 
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of Eq. (^|) at a much later time than is considered in Eq. (|10D . Consider the ensemble average 
of Eq.(|9p, but now at times for which the nuclear fields are time dependent. In addition, 
we average this quantity over a time interval large compared to the period of the electron 
precession but small compared to the time at which the ensemble average is taken. This 
leads to: 

(g(t» = (n {t )(n(t) ■ S(t))) = ( B " mB B ^ ' m) ) ■ (12) 

In this equation (Bjf(t) ■ S(t)) is the energy of the electron-nucleus spin system, which does 
not depend on time and is equal to (B N (0) ■ S(0)). Nor does 7?^ change its value, because 
dB^/dt _L Bjq. As a result, 



m)= ( Myi ) . (13) 



The rotational symmetry of the Hamiltonian Eq.|] leads to (BN{t) a BN{0)f3) = 
6 a> p(B N {t) x B N {Q) x ) = {6 atP /3)({B N (t) -B N {0))). This allows us to write Eq.|13| as 

We see again that (S(t)) is directed along S(0), which is the only physically defined direction 
in the ensemble. This equation describes the time dependence of the electron spin relaxation 
at times large compared to 7a and is limited only by the characteristic time scale of the 
dipole-dipole interactions. One can see that at times < t < T N Eq.(|TJ]) gives the same 
result as Eq.(|T0|). 

The variation of nuclear magnetic field direction in the time dependent correlation, in 
Eq. (|HD , is limited by the conservation of the total spin angular momentum of the electron- 
nucleus spin system , F = S+J2j I ■ The conservation of F at times shorter than the nuclear 
dipole-dipole relaxation time follows from the fact that F commutes with the Hamiltonian of 
t e electron-nucleus spin system Eq.(E3). The total nuclear spin is, then, effectively conserved 



since: 7s = Z)j Ij = F — S ~ F because F ^> S The latter follows from the fact that 
the dispersion of 7^ increases with the number of nuclei in each QD and the average value 
of |7 S | > 1. 
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If the electron wave function were constant in the localization region and zero outside (the 



so called "box model" [jHJ) the nuclear magnetic field, B N , would be proportional to the total 
nuclei spin, Jg, and would also be conserved. From Eq.(|T4|) the nuclear spin precession would 
not then lead to any additional spin relaxation in this model, and (S(t)) = So/3- However, 
in real QDs the amplitude of the electron wave function at the nuclei in the localization 
region does depend on their position. The nuclear field is not uniquely determined by the 
value and direction of the total nuclear spin; a distribution of values of Bn are possible for 
the same value of is- The total nuclear spin, Js, can be distributed in different ways among 
the nuclei as a result of their interaction with the electron. 

The determination of the correlation, Eq.(|14D at all times is beyond the scope of this 
paper. We will evaluate it only in the limit t ^> Tjq, where the B^{t) are randomly dis- 
tributed. The joint distribution function of the nuclear fields and total nuclear spin can be 
written as the product of the total nuclear spin distribution function and the conditional 
probability distribution of nuclear fields given a certain value of the total nuclear spin: 



W{B n ,Iy) = W i {Iy)w(B n \I t ) . 



(15) 



Both of these latter distributions have a Gaussian form [[U; 

1 



7r3/2 A 3 



cxp 



:ir) 2 



A? 



(16) 



and 



w(B N \Iz) = , , 3 exp 

7T / Ajj 



[B 



N 



A' 2 



(17) 



where A 2 - = 2NI(I + l)/3 is the dispersion of the total nuclear spin distribution, and iV 
is the number of nuclei which are in statistical equilibrium with the localized electron (see 
Appendix). This number is approximately equal to the number of nuclei in the dot. The 
conditional distribution of the random hyperfine fields, given a value is for the total nuclear 
spin, is shifted from zero to the field due to the weighed total nuclei spin, (a) is, where 



a) = (£f dj)/N = A/(n B g e N). The dispersion of this distribution, AJ = A| - (a) 2 Af, is 



determined by integrating Eq.[l7|with distribution Eq.(|l6D over all Js and comparing it with 
Eq.(|[). Clearly A J < A 2 B , because the latter also includes the dispersion of the random 
distribution of the total nuclear spin. 

In general, to calculate the correlation between the direction of the hyperfine nuclear 
field at two times, n = B N (0)/B N (0), and n(t) = B N (t)/B N (t), (see Eq.flPD), we have 
to find the angular distribution of the magnetic field Bjv(t) given a certain value both of 
the total nuclear spin J and the magnetic field magnitude, B^ = |-Bjv|- This conditional 
distribution function is written: 



F(n\B N ,I) 



exp 


-(B N n-(a)I) 2 /A'i 




Jdn(n) 


exp 


'-(B N n-(a)I) 2 /A' 2 



(18) 



Using distributions of Eqs. (^, p!6| , |18D we find the correlation Eq.(|TJ) 

7 = ((n-n )) 

n ■ n )F(n\B Nl I)w(B N \I)W I (I)d 3 B N d^l(n)d 3 I . 



(19) 



This gives us the long time (t 3> T/v) electron spin polarization (Eq. (|14D . Straightforward 
calculation gives: 



irx 3 Jo 



dy 



dz 



[2y cosh(2y) - sinh(27/)]^ 



x exp 



1 + — )z 2 - - 
X 2 J \z 



yz sinh(2y) 

2" 



(20) 



where 



x 



(afA^ 
A'j 



(a> s 



(21) 



(a 2 ) - 

is the relative dispersion of the hyperfine magnetic field of the localized electron acting on 
nuclei and (a 2 ) = Ajj/A 2 -. If all the nuclei in the unit cell have the same hyperfine constant 
x = N L /(N-N L ). 

In Fig. 1. we show 7 as function of (a 2 ) /(a) 2 ~ Nl/N. In the limit of large x, where 
the number of nuclei contributing to the field, Bn, is close to the number of nuclei whose 
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total nuclear spin is conserved, j(x) ~ 1. This case is realized in the "box model" and also 
at short times t < Tjv- The number of nuclei that are included in the electron-nucleus spin 
system of the Hamiltonian Eq.(2), and whose total spin is conserved, increases with time, 
which in turn decreases 7. However this only holds for t < Td-d the dipole-dipole nuclear 
spin relaxation time in the dipole field of its neighbors. For spherical QDs with 10 5 nuclei 
confined in an infinite potential barrier (note that Nl ~ N tot /2.8, where A^ot is the total 
number of nuclei in the QD), our approach is valid for the time interval 10T/v < t < 100T/v 
(see Appendix). The calculation in this model shows that the number of nuclei that need be 
included at t ~ Td-d is ~ 2.1A/, (7 ~ 0.4) and does not reach its maximum value N « 2.8Nl. 

It is important to note that non vanishing average spin polarization at time t ^> Tjv 
(7 7^ 0) means that there is a significant probability that an ensemble of electron-nucleus 
spin systems retains its initial spin state. 

IV. SPIN DECOHERENCE IN A STRONG EXTERNAL MAGNETIC FIELD. 

A strong external magnetic field, B, (B 3> B^) significantly changes the process of 
electron spin relaxation. In this large field the Zeeman splitting of the electron spin levels 
is larger than their inhomogeneous broadening in the hyperfine nuclear magnetic field. The 
total magnetic field acting on the electron is now effectively directed along the external 
magnetic field. The nuclear hyperfine fields only perturb the precession frequency of the 
electron spin about the external magnetic field direction. 

Consider, now, the effect of a strong external magnetic field on the electron spin polar- 
ization. The motion of the spin in the total magnetic field is again described by Eq. (|9|) 
where, now, n = (B + B N )/\B + B N \. Averaging Eq. (0) over the ensemble, using the 
distribution of nuclear magnetic fields in Eq.(^), we obtain: 

(S(t)) = R\\(t)(S ■ b)b + R ± (t)[S - (S ■ b)b] 

+ R 1 ± (t)[(S -(S -b)b)xb] , (22) 

where 

11 



i?n(t) = + Ai2||(t) , 
= #j° + A#}_(t) , 

and b = B/B is a unit vector along the external magnetic field, 

/ [S + (B Jv .b)] 2 \ 



and 



iff 



1 / B%-(B N -b) 2 ' 
2\ |S + BrJ 2 



(23) 



— f°° ,1 [°° j v ex p(- z2 -y) 

sR Loo Jo y (B + zY + y 



(24) 



2y/7rJ-oo Jo (3 + z)* + y 
is the value of (S(i)) in the long time limit £ ^> T A , and /3 = B/A B . The time dependent 
components are given by 



AR\\(t) 



OO />00 



dz I dy 



yexp(-z 2 -y) 



t/tt J-oo Jo ((3 + z) 2 +y 
X cos 



AR° ± (t) = E(t) - 



y/(P + z) 2 + y 
AR l{ (t) 



t 



(25) 



where 



and 



E(t) = exp 



( t 



V2Ta 



f(3t\ t . (f3t\ 

cos — — — — sin — — 

\T A 2(3T A \T A 



p i m 1 r°° , r A qi M + z ^ ± z vY^- z2 - y) 
H±{t) ~V^J-oo Jo y WTzfTW 2 



x sin 



(P + zy + y^- 

1 A 



(26) 



(27) 



Equation ( P2"D simplifies considerably for strong magnetic fields. Calculation of the coeffi- 
cients in Eqs. (|23|j24lj25|) in the limit /3 3> 1 gives: 



(S(f)> « 

f 1 - cos (u B t) 
1 ^ exp 



\2T A J 



:Sn • b)b 
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+ 



. . 1 — cos (uJ B t) 
cos{u B t) + — 2 



+ sin(^ B t)[(5 - (Sq ■ 6)6) x 6]}exp 



[So - (So ■ 6)6] 

t 



(28) 



,2T A . 

One can see that in strong magnetic fields, B ^> A# the component of spin along B is 
conserved, while its two transverse components precess with a frequency ujb = I^bQcB /h, 
and decay as a result of the inhomogeneous broadening of the levels in the random magnetic 
field of nuclei, respectively . The dephasing arises from the dispersion of the nuclear field 
along the external magnetic field, which leads to an inhomogeneous dispersion of the electron 
precession frequency. The perpendicular components of the nuclear magnetic field change 
the direction of the precession axis by a small angle of ~ B^jB and lead to a dephasing 
rate that is (B^/B) 2 smaller than that due to the dispersion of the nuclear field along the 
external field. This result is consistent with the small ratio of the lifetime broadening to 
the secular broadening describing the transverse relaxation time T2 in strong magnetic fields 
(see e.g. Ref. |7|). 

Figure 2 shows the time dependence of the various components of (S(t)) which occur in 
Eq.(^2]). The both longitudinal and transverse components of the electron spin polarization 
tend to a steady state value after several oscillations. The number of oscillations grows with 
increasing a magnetic field. Increasing the magnetic field also changes the steady state value 
of the longitudinal component i2ii(oo) (from 1/3 to 1) and the transverse component R° ± (oo) 
(from 1/3 to 0), while the steady state value of R^oo) is zero for all values of B. In a 
strong magnetic field when the nuclear spin relaxation mechanism of the longitudinal spin 



polarization is suppressed phonons can again play an important role (see ref. ||18|| ). 

The important characteristic measured in steady state experiments, such as Hanle 
effect measurements (e.g., see ||), is the average electron polarization for its lifetime, 
(1/t) J(S(t)) exp(—t/r)dt, where r is the lifetime of the localized electrons. Comparison 
with Eq.(|22|) shows that this average polarization is characterized by: 



1 r°° 

P][ (P,t) = - R\\{t)exp{-t/T)dt , 
r jo 
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1 r°° 

p°/(P, t) = - R°/(t) eM-t/r)dt . (29) 
r Jo 

The dependence of these respective terms on the magnetic field is shown in Fig. 3 for 
different values of the electron lifetime. 

Equation (122|) , as is Eq. ([To|) , is derived assuming a time independent frozen fluctua- 
tion of the nuclear field acting on the electrons. In zero external magnetic field, this no 
longer holds at longer times such that we must take into account the nuclear spin pre- 
cession in the inhomogeneous hyperfine field of the electron. Nuclear spins precessing at 
different rates about the average electron spin direction create a time dependent hyper- 
fine field with components perpendicular to the original direction of the frozen fluctuation. 
As a result the average electron spin projection follows the new direction of the slowly 
varying nuclei field. The characteristic time of this slow process is determined by the dis- 
persion of the nuclear spin precession frequency in the inhomogeneous field of the electron 



Ttf 1 ~ (n B ge/K)\l (a?) — (a) 2 ~ {^Bge/K)\J {a?) . This second regime of spin relaxation be- 
gins when t > Tjy. 

In a strong external magnetic field B ^> A#, the average electron spin is directed along 
this strong field, independent of the nuclear hyperfine fields (Bn <C B). Although the 
nuclei precess with different frequencies in the inhomogeneous electron field, the electron 
is effected only by the component of the nuclear field along the external field. As a result 
the nuclear magnetic field acting on the electron spin is frozen for times much longer than 
Tjy. Thus a frozen fluctuation model of the nuclear hyperfine field is valid when describing 
the dephasing dynamics of the electron spin polarization in an ensemble of quantum dots in 
strong magnetic fields. As we mentioned above this consideration is limited by a low enough 
temperature and the time scale of the nuclear dipole-dipole interaction. 

In each QD, the motion of the electron spin in the hyperfine field of the frozen fluctuation 
of the nuclei is coherent. The dephasing is a result of inhomogeneous broadening of the 
electron spin levels in the ensemble of quantum dots. This makes it possible to recover the 
transverse electron spin polarization using the spin echo technique |T^|, which also can be 



14 



used for quantum computation [|T^ . 



V. DISCUSSION 

We have determined the time dependence of the electron spin relaxation rate arising 
from its interaction with nuclear spins for an ensemble of QDs, or equivalently, averaged 
over a large number of successive measurements of a single dot. This gives us the value of 
T 2 * usually measured in ensemble experiments. 

The electron spin relaxation time was studied in three types of QDs: GaAs natural 
quantum dots, which are island-like dots formed by a fluctuation of the GaAs quantum 



well thickness ||20|| ; self organized InAs QDs pl| , and nanocrystal CdSe QDs [P2j| . We can 
estimate the value of for each of these dots. In GaAs natural QDs, using values of A 7 
from Ref. [[13] and Eq.(6), and taking the dot to contain 10 5 nuclei, Eq.(ll) gives Ta ~ 1.0ns 
The values of A- 7 are not experimentally determined for InAs. Assuming that compound 
ionicity does not significantly vary among the semiconductors GaAs, InSb and InAs ||T3| , |2"3"|| , 
we take the hyperfine constants, A- 7 , for As and In ions from |13|] and ||24|| : A 3 As = 47/j,eV and 
A J In = 56/ieV, respectively. The large value of the In nuclear spin, I = 9/2, strongly effects 
the magnitude of the hyperfine interaction in InAs, and shortens the spin dephasing time. 
To estimate this time we need to know the number of nuclei in the electron localization 
volume. In the self organized QDs, where the electron wave function strongly depends on 
the QD shape and barrier height, we used the geometrical volume of the dot to determine 



the number of nuclei, Nl, in Eq.(7). The QDs studied in |Tl|] contain 1000-4000 nuclei 



25] . This leads to electron spin relaxation times on the order of 50-100 ps, which is close 
to the value of gT 2 * measured in these dots. The hyperfine interaction constants in CdSe 
are also not experimentally determined. In this material only a fraction of the nuclei (25 
% of the Cd ions) have a magnetic moment and these have spin 1 = 1/2. As a result, 
the electron spin interaction with the nuclei in these QDs is weak. We shall take for the 
hyperfine constant in CdSe half that of Cd 111 PBj, assuming completely covalent bonding 
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in this material. This leads to A 3 Cd = 12/^eV. Taking the CdSe lattice constant a® pa 4.2A 



[27 1 we calculate Nl = 3, 500 for the ground IS electron state in spherical nanocrystals with 
a radius 28A. Only 1/8 of these nuclei contribute to dephasing. This leads to Ta ~ 1.6 ns, 
which is consistent with the experimental measurement T 2 * = 2.5 ns [^2| 

In conclusion, we developed a theory of electron spin relaxation in QDs arising from 
their hyperfine interaction with the QD nuclei. The relaxation is determined by three 
physical processes: (i) - the precession of the electron spin in the hyperfine field of the 
frozen fluctuation of the nuclear spins; (ii) - the precession of the nuclear spins in the 
hyperfine field of the electron; and (iii) - the precession of a nuclear spin in the dipole 
field of its neighbors. These processes have three disparate characteristic times. For GaAs 
QDs with 10 5 nuclei, they are ~ Ins, ~ 300ns, and ~ 100/is, respectively. The last of 
these times is so long that many other electron spin relaxation mechanisms can be more 
important on this time scale. An external magnetic field suppresses the relaxation of the 
spin component along the magnetic field. The transverse components of the electron spin 
polarization relax completely in a time on the order of the electron precession period in the 
field of the frozen nuclear fluctuation. Comparison with experimental data shows that the 
hyperfine interaction with nuclei is the dominant mechanism of electron spin relaxation in 
QDs. 
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APPENDIX: THE NUMBER OF NUCLEI INTERACTING WITH A LOCALIZED 

ELECTRON IN A QD 

The number of nuclei, N, "interacting" with a localized electron and forming an electron- 
nuclear system in thermodynamic equilibrium is generally a time dependent quantity. For- 
mally, the electron interacts at all times with all those nuclei of the crystal lattice sites where 
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its wave function is nonzero. However, the strength of this interaction decreases with the 
distance from the QD center as the square of the electron wave function, and the random 
magnetic field Bjv, which determines the electron spin precession is determined only by the 
relatively small number of nuclei, N L) in the" effective" volume of the electron localization 
V L (see Eq.8). This volume (or alternatively N L ) completely determines the electron spin re- 
laxation at short times t < Tn during which nuclear spins, and, therefore, the corresponding 
field B N do not change their directions. 

The precession of the nuclear spins in the hyperfine field of the electron causes time 
dependent changes in B^. This time dependent B^{t) in each dot is determined by a 
macroscopic number of initial conditions. The correlation (B N (t)B N (0)} (see Eq.14) is now 
an important characteristic of an ensemble of quantum dots. The ensemble of localized 
electrons completely loses its memory of the initial spin polarization if this correlation is 
zero. Although (B N (t)B N (0)) can be calculated for any time, t, in principle, in the long 
time limit t ^> T N , when the values of B N (t) are randomly distributed, one can find this 
correlation using statistical physics techniques. 

Let us separate the nuclei interacting with a localized electron during a time t into two 
groups: (1) those nuclei whose spin rotates through an angle cutf > 1 and (2) those whose 
spin rotates through an angle uj{t < 1, where Ui is the precession frequency of the nuclear 
spin in the hyperfine field of the localized electron. The second group can be considered 
as nuclei that "do not interact" with the electron. For times t 3> T N the first group 
automatically includes the Nl nuclei that determines the field B^. The number of nuclei 
N(t) "interacting" with the electron, and thus belong to the first group, increases with time 
and determines B^it), whose variation is restricted by conservation of the total spin of 
the electron-nuclear spin system. Randomness of the distribution of the N(t) nuclear spins 
limited by the conservation of B N and J s is the one major assumption of this statistical 
model. 

For nuclei that interact strongly through the "indirect" hyperfine interaction, i.e. nuclei 
whose spins rotate through an angle uit ^> 1 in the electron hyperfine field, this assumption 
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is justified. The error in the statistical approach is connected with the border nuclei for 
which u{t ~ 1. These nuclei strongly interact with the electron, but, on the other hand, 
they still "remember" their initial direction at t — 0. The number of such nuclei is on the 
order of AiV pa tdN / dt. As a result the relative error of our estimate is SN = AN(t) /N(t) ~ 
(t/N)dN/dt. 

The time dependence of N(t)/NL and SN in spherical QDs with an infinite barrier is 
shown in Fig.4. The time is measured in the unit of TjJ- ~ (32 / V) (h/ A) (a 3 /vo), the nuclear 
precession time at the center of the QD, where a is the QD radius. In GaAs QDs with 10 5 
nuclei, this time is on the order of microseconds. One sees that N < Nl for t < lOTjy. In 
this time period 5N > 1 and our statistical approach is not valid. However N(t) increases 
and 5N(t) decreases rapidly with time, and SN ~ 0.25 at t pa 20T^. After this time the 
nuclei which give the main contribution to the magnetic field, B^, acting on the electron 
are in statistical equilibrium. This allows us to use Eq.(|20"D. 
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FIGURES 

FIG. 1. The dependence of the ensemble averaged electron spin polarization on the ratio of 
dispersion of the electron hyperfine field acting on the nuclei to its ensemble averaged value. 

FIG. 2. The time dependence of the longitudinal, (a), and transverse, (b) and (c), components 
of the ensemble averaged electron spin polarization calculated for different magnetic fields. The 
curves are calculated for (3 = 0, 1, 2, 3, 5 respectively. 

FIG. 3. The magnetic field dependence of the longitudinal, (a), and transverse, (6) and (c), 
components of the steady state electron spin polarization. Calculations are done for electron 
lifetimes r/T A = 0, 0.2, 0.4, 0.8, 1.4, 3.5, and 10. 

FIG. 4. Time dependence of the ratio of the number of nuclei, N(t), "interacting" with an 
electron, in a spherical QD with an infinite potential barrier, to Nl and the relative error, SN(t). 
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